in the same limit in which RPA itself is exact. 
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Based on an exactly solvable model of a nearly homogeneous solid (jellium+phonons) , 
we compute the optical dielectric function exactly at zero temperature in any number of 
spatial dimensions. It is argued as we have done elsewhere, that the answer for the optical 
■ dielectric function we obtain is in fact the exact answer for a real nearly homogeneous solid 
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I. INTRODUCTION 

In earlier preprints0, we introduced a scheme that involved replacing the electrons in a solid with sea- 
bosons. This has the distinct advantage of rendering any problem that involves interacting electrons, exactly 
soluble in the RPA limit. In this article, we use the same approach to compute exactly the optical dielectric 
t^"- ■ function of a nearly homogeneous solid. 
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II. MODEL HAMILTONIAN 

a 

I 1 The model hamiltonian of a nearly homogeneous solid may be written as, 

H = 5> k (q)4(qH(q) + £ ^ £[A k (q)a k (-q) + A k (-q)a k (q)][A k , (-q)a k , (q) + A k , (q)a£, (-q)] 

^ ' k,q q^O k k ' 
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Here as usual, 



and, 



E 77T £[A k (qW(-q) + A k (-q)4(q)](6 q + bL q ) + ]T fi q & q 6 q (1) 



q 



Wk(q) = A k (-q)^ (2) 



in 



A k (-q) = v /n F (k-q/2)(l-n F (k + q/2)) (3) 

and, np(\<L) — 0(kt — |k|). This hamiltonian consists of the usual coulomb terms and the new ingredient here 
is the interaction with phonons. We assume that there is only one branch for simplicity, the optical phonon 
branch. In order to diagonalise this, we proceed as usual, by introducing a Bogoliubov transformation^]. 
Assume that the diagonalised form looks as follows, 
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q,i 

The relation between the dressed sea-displacements d qi , and the parent bosons may be deduced as follows, 
The dressed sea-displacements and the original sea-displacements may be related by, 

d ^ = E^' 4(q)] a k(q) - E^' a k(-q)] ft k(-q) + b q ]& q - &- q ]&t q (5) 

k k 

Ok(q) = X)K(q)> 4q,iK,i - Xt ak ( q )' d -q.il d - q ,i ( 6 ) 

&q = EK dq.Jdq.i ~ XK d-q.d^-q.i ( 7 ) 

K,i,a k (-q)] = -( - ^ k(q) )(%i)G 2 (q,.) (8) 



here. 



K,-4(q)l = ( Ak( q ) n (% i )g2(q,0 0) 



K,^ q ] = (^)( ^-)G 2 (q,i) (10) 



K,i^-q] = -(^)(— ^ r )G 2 (q,z) (11) 

w q,i + "q 



GM = [(^) ( ^^| )2 +g -^( q ,0(y) 2 ]-* (12) 
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gq,i=« q + Mg „ ^ (13) 
w q,» "q 



5 - 2 (q, = Et ^^'^Tq!^^ ] (14) 
k m > 

The eigenvalue w q! , is given by, 

£i?PA(q, w) = o (15) 

and the RPA dielectric function with phonons is given by(analytic continuation is by replacing i uj — > 
w + i 0+). 

i ■ \ , ,-5 q . w , n F (k + q/2) -n F (k-q/2) 

e fl P A (q,*w) = l + H*r)E ■ — (16) 
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K-^hyl (17) 



The optical dielectric function is given by, 

e(u>) = 1 - =| - ^7r(w) ( I S ) 

where 7r(w) is the retarded form of the current-current correlation function^ : 
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?r(a;) = --/ dtQ(t-t) exp(iuj (t - 1 ))([3(t),3(t )]) (19) 
The total current is given by, 

J = e V A k (-q)(-i) 4(q)a k (q) (20) 

* — ' TO 
k,q 

If the current operator J is independent of time then tt(uj) = 0, and the imaginary part of the optical dielectric 
function is identically zero, indicating that there is no absorption of radiation. This is going to happen in 
the case when phonons are absent and then we know that the total current is conserved in a homogeneous 
electron-gas. This fact is obvious in the case when we use the undiagonalised form, but less so when we use 
the dressed-sea displacements. In order to verify this latter fact(namely that J is time independent in the 
case when phonons are absent), let us write down the relevant quantity (possibly time dependent part of the 
total current) and argue that it is zero, 

Part of J = -e XI H[ ak ( q )' K(q), d -^M\^-^3 exp ^ ( w q,* + w q,j)*) ( 21 ) 

q k i,j 

This is going to be zero if there is only one root to the eigenvalue equation, Eq.( [l5]). This is the case when 
there are no phonons present. Even in more than one dimension, without phonons, the equation has only 
one root. This is the alternative (somewhat roundabout) explanation as to why the homogeneous electron 
gas absorbs no radiation. But when there are phonons present, the eigenvalue equation has precisely two 
roots (no matter what the dimension is). This is why J is now time dependent and this leads to a finite 



imaginary part of 7r(w). In order to demonstrate that there are precisely two roots of the equation Eq.( 15), 
it is sufficient to solve the ID case as this is likely to be indicative of the state of affairs in more than one 
dimension as well. Unfortunately, the solution is not simple. One has to solve recursively the equation, 



* (— h 

to 



The (postive) imaginary part of this root is the eigenvalue we are looking for. The recursive solution obtained 
by replacing v qjU sucessively by first v q , and then later on by v qMo where loq is the zeroth order root and so 
on, is one such root. The other root is obtained by recursively solving the same equation differently. Let us 
rewrite the recursion as, 
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/ 2 M2fi„ 
u = ijci 2 „ 5 q — 23 

where, 

S( ,, w) _l±<^zs)!V ,24, 

This may also be solved recursively, the zeroth order approximation being, iQ q and the rest as discussed 
before. Thus we have eigenvalues, in the weak coupling regime they are close to kf\q\/m and £l q respectively. 
But in the strong coupling regime the answers are going to be different. The situation is likely to be similar 
in more than one dimension. Let us now compute the time evolved total current. 

J(t) = e £(-)M-q)4 *(q)4(q) (25) 

' TO 



k.q 



now, 



4(q) = E MqMq.Jdq 

i=l,2 

- X>k(q),^q,<]d t _ qii e iw *' t (26) 



Therefore, 



([J(t),J(t )]) = ^ F(q ;l ,j) S m[K,, + Wqj )(t' -i)] (27) 
^(q;i,j) = (e 2 )(2 0(^)X>k(-q)A^(-q)[a k (q),d_ q Ja^ 

k,k' 

- (e 2 )(2 i)(^) £ A k(-q)A k ' (q)K(q), d- q> i][a k (q),< jJ ][a k ' (-q), ^ q JK' HO. dqj] (28) 

k ,k' 

This leads us to the following formula for the retarded form of the current-current correlation function, 

««) = E E(%^) ^— -7F (») 

q i^j 

The formula that relates the refractive index and the absorption coefficient to the optical dielectric function 
has been given in the review by Haug and Schmitt-Rinki. 

/ x i , s c a(uj) 
e(w)s = n(w) + i _ v ; (30) 

The imaginary and real parts are of e(u;) are easy to write down, 



q>2 ^qj 



The refractive index is given by, 

,1 



,(,:) : : [(~)[Re(e(u;)) + {{Re(e{u))f + (Zm(e( W ))) a )>]]> (33) 

q; |/m(e(a>))| 
n[cj) c 
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